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Abstract. We discuss the experimental feasibility of quantum simula- 
tion with trapped ion crystals, using magnetic field gradients. We de- 
scribe a micro structured planar ion trap, which contains a central wire 
loop generating a strong magnetic gradient of about 20 T/m in an ion 
crystal held about 160 /im above the surface. On the theoretical side, 
we extend a proposal about spin-spin interactions via magnetic gradi- 
ent induced coupling (MAGIC) [Johanning, et al, J. Phys. B: At. Mol. 
Opt. Phys. 42, (2009) 154009]. We describe aspects where planar ion 
traps promise novel physics: Spin-spin coupling strengths of transversal 
eigenmodes exhibit significant advantages over the coupling schemes in 
longitudinal direction that have been previously investigated. With a 
chip device and a magnetic field coil with small inductance, a resonant 
enhancement of magnetic spin forces through the application of alter- 
nating magnetic field gradients is proposed. Such resonantly enhanced 
spin-spin coupling may be used, for instance, to create Schrodinger cat 
states. Finally we investigate magnetic gradient interactions in two- 
dimensional ion crystals, and discuss frustration effects in such two- 
dimensional arrangements. 



Introduction 

When many atoms of equal or different species bind together and form a piece of 
matter, its properties depend on the mutual interactions of the constituents and the 
structure which was formed. Surprising collective phenomena appear in the magnetic 
properties of a lattice, for example when spins are aligned parallel or anti-parallel with 
each other, as observed in ferromagnetic or anti- ferromagnetic materials, respectively. 
Also the electric transport properties in such materials can undergo phase transitions, 
or insulators or superconductors may form, depending on control parameters such as 
chemical composition, dopant concentrations, pressure, temperature and magnetic 
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field. The goal of predicting such material properties from first principles and even 
to design these properties in view of applications is difficult because the complexity 
of quantum interactions increases rapidly with the number of particles. Here, the 
importance of quantum optical model systems becomes apparent, as one can construct 
a complex system from individual, fully controlled atoms or ions, where even the 
mutual interactions are tunable by external control parameters over a large range. 
This bottom-up approach of simple model systems promises to allow a verification 
of the computational methods, as well as the investigation of the basic properties of 
quantum magnetism and superconductivity. A long term goal is thus the simulation 
of real material systems with such quantum optical model systems, following the 
vision of R. Feynman who, in 1982, proposed simulating the static and dynamical 
properties of systems that are difficult or even impossible to solve numerically by 
investigating experimentally a well controlled model system |1J. This upcoming field 
has been attractive for both experimental and theoretical physicists |2|3j . 

Experiments with trapped atoms in optical lattices have pioneered the field, follow- 
ing the proposal by Jaksch, Cirac and Zoller [4.,5,6j with spectacular demonstrations of 
the Mott insulator to superconducting phase transition [7] . Today these experiments 
aim for complex lattice structures, including superlattices, frustration in Kagome 
arrangements [8 9 , or a simulation of SU(N) magnetism [lQ|llj. Trapped ion quan- 
tum simulation has been proposed by Porras and Cirac [12 , where optical dipole 
forces from laser beams induce spin dependent forces to generate effective spin-spin 
interactions. Experimentally, we have seen the simulation of a magnetic quantum 
phase transition with two ions [13 and more recently the simulation of a frustrated 
three-spin anti- ferromagnetic interaction with three ions [14 • Designed Hamiltonians 
may be used to simulate the Dirac equation and relativistic Hamiltonians |15|16j , or 
crystal-to-superfluid transitions in XXZ-spin chains [TT. Yet another way to design 
spin-spin interactions has been proposed by Wunderlich [18 , where strong static mag- 
netic gradients are used. In an inhomogeneous magnetic field, spins experience forces 
such that the ions are displaced from their equilibrium positions in the ion crystal. 
This leads to a mutual effective spin-spin interaction. First experiments with a single 
ion have demonstrated the spin-dependent force [19 ; the generation of many-particle 
entangled quantum states [20 and a generalization with rf-magnetic fields [21|22j 
has also been proposed. As the interactions are generated without laser light fields, 
but only by magnetic gradients, scaling of the simulation to more particles might be 
facilitated [23 . 

We describe in this paper how magnetic interactions may be simulated when a 
cold crystal of ions in a planar ion trap is used. We discuss the specific advantages of 
our approach in a particular two-dimensional trap, where ions come close to the sur- 
face, whereby they may experience large magnetic gradients which are necessary for 
significant spin-spin interaction strengths. Furthermore, the planar trapping device 
should permit the arrangement of ions for a quantum simulator either with ions in one 
linear crystal or in multiple parallel linear crystals, thus forming a two-dimensional 
structure. In Sec. |1.1| we will therefore summarize the basics of planar traps and de- 
scribe the specific features of the device we have developed. We estimate the strength 
of the expected magnetic field gradients in Sec. |1.3| Sec. [2] then shows how spin-spin 
interactions are generated and how the strength can be designed depending on the 
applied magnetic field gradients and their direction, allowing one to switch from a 
ferro- to antiferromagnetic interaction. Sect. 2.4 describes how the time dependence 
of the corresponding current resonantly selects specific eigenmodes of the crystal. 
A scheme for a two dimensional ion crystal well suited for quantum simulations is 
presented in Sec. [3j as well as the complete calculation of eigenmodes and spin-spin 
interaction strengths in Sec. 3.2 All parameters given in the paper are within reach 
of our on-going experiment using the particular planar trap [24 also described in 
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Fig. 1. Typical arrangement of the electrodes 
in a planar trap. A linear trapped ion crystal 
is sketched oriented along the z-axis. 



Sect. |l.l[ The goal of this paper is to highlight possible avenues towards the quantum 
simulation of magnetic systems with ion crystals and to stimulate the discussion of 
further types of spin interacting multi-particle systems. 



1 Planar traps with magnetic field gradients 
1.1 Planar micro traps 

Using static and dynamic electric fields to hold ions, the planar trap evolved from the 
linear Paul trap [25 . Contrary to the latter, where the electrodes are configured in 
three dimensions, in a planar trap they are all aligned in one plane and the ions are 
trapped above the surface housing the electrodes. This adaptation from its previous 
form allows easier microfabrication using advanced lithographic techniques, which is 
especially important when aiming for the complex structures required for scalability 
towards large numbers of ions [26]. On the other hand, in planar traps, where the ions 
are typically 40 — 200 jam above the conducting electrodes, the potential is shallow 
and the close vicinity to the metallic surface may lead to heating |27|28|29l . Most of 
the up-to-date realizations utilize five parallel electrodes [30'31' (see Fig. [l|). Two of 
them are driven by radio frequency (rf) voltages creating a quadrupolar dynamical 
trapping field in x and while the others are held at static voltages (dc) which 
confine ions in the z direction. 

The resulting potential for a particle with charge q and mass m is t) = 

^dc(^7 ^, z) + 0rf z) cos Qt leading to an equation of motion given by 



(Pxi 



9>dc ^ , f d^rf 



dxidxj I \dxidx 



cos{f2t) 



0, (1) 



where f2 is the angular frequency of the rf potential, xi = x, X2 = xs = z and the 
potential is expanded to second order around the equilibrium position. Rewritten in 
standard Matthieu form, the equation of motion is 

mj-^ + [A + 2Q cos (2r)] x = 0, (2) 



where 



mfi^ \dxidxj J ' mfi^ \dxidxj J ' ^ ' 

are the entries for matrices A (stiffness matrix) and Q (excitation matrix) determining 
the stability parameters of the trap [32 and r = Qt/2. The trajectory solution of 
the ion undergoes oscillations including both the so-called micromotion frequency at 
Q /2ti and the secular motion frequencies at where ~ \/ An + Q%l2 in the 

case that the matrices are diagonal and the entries are much smaller than one. For 
the accurate numerical calculation of the trap potential without approximations we 



4 



Will be inserted by the editor 




Fig. 2. Left: Top view of the planar trap with optical microscopy ^24^. The 5 wire configu- 
ration with dc, rf, dc, rf, dc, is modified with the center dc electrode divided into two wires 
Wi and W2 and a nine- fold segmentation of both outer dc electrode. The whole chip has a 
width of about 10 mm. Right: Zoom of the center of the trap. The outer dc electrodes have 
a width of 595 /xm, the rf electrodes of 197 /xm and the inner dc electrode Wi a width of 100 
/xm. The gold layer is 1.9 /xm thick. The gap between rf electrodes and their neighbors is 
12 /xm to avoid electrical breakdown when applying high voltages. A gap of 3 /xm is chosen 
between dc electrodes. 




Fig. 3. Scanning electron microscope pic- 
ture showing the quality of the insulating 
gap between the gold electrodes. The scale 
bar is 1 /xm long. The surface roughness has 
an RMS value of less than 5 nm, and a peak- 
to-peak edge roughness of 40 nm. The gold 
resistivity is 2.63 /xi7cm. 



employ a numerical toolbox [33] reaching a 1 % accuracy [34]. With the geometry of 
the present chip (see Fig. [2]) we calculate that the trap center is yo = 164 /xm above the 
surface. Along the z-axis the segmented dc electrodes allow trapping both above the 
center of the winding inner electrode and also at a distance of more than ± 1000 /xm. 
Using an rf amplitude of Vi-f = 230 V at a frequency of i? = 27r • 34.88 MHz the depth 
of the trapping potential for ^^Ca^ is 150 meV and the trapping frequencies are Ux = 
27T • 3.14 MHz and ujy = 27r ■ 3.20 MHz in the radial directions, and cj^ = 27r • 310 kHz 
in the axial direction (see Tab. 1 in Appendix B). 

For the fabrication of such devices, various methods may be used. Photolithogra- 
phy, electroplating and etching on Borosilicate glass has led to a first generation of 
traps. For the details of the fabrication process we refer to Ref. [35]. A newly estab- 
lished process builds on the good heat conductivity and surface quality of a sapphire 
substrate, on which a photolithographic procedure is applied. The developed resist 
serves as a mask for the deposition of a 50 nm thick adhesion layer of titanium and 
a high quality gold film of about 2 /xm thickness, both being deposited by a ther- 
mal evaporator. The Ti and Au deposited on the resist top are removed in a lift-off 
process. From an AFM analysis we find an rms surface roughness of less than 5 nm, 
while an SEM analysis determines the edge roughness to be 40 nm (peak-to-peak) 
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Fig. 4. Trapped linear ion crystals of Ca^. The applied radio frequency at 30.82 MHz with a 
Vpp=400 V and the dc voltages of 2.2, 0.7, -0.8, 0.7, 2.2 V on five adjacent pairs of electrodes 
result in an axial frequencies of uoz — 27r 234 kHz. 

(see Fig.[3|. The process has an overall alignment accuracy of better than 1 /im. In 
our tests, we continuously apply 205 V amplitude without observing any sparking 
which we attribute to the high quality of the electrodes without sharp edges. 

1.2 Trapped linear ion crystals 

The above structure has been loaded with linear ion crystals of ^^Ca+, see Fig. |4j 
We load the ions by a photoionization with laser light near 423 nm and 375 nm. 
The ions are cooled on the dipole transitions near 397 nm and 866 nm. The emitted 
fluorescence is imaged on a EMCCD camera. 

1.3 Integrated magnetic field gradient coils 

The trap design includes a split inner electrode with wires Wi and W2 (see Fig.|2|. 
The small loop in Wi with inner length 150 /im and outer length 250 /im generates 
a high gradient of the magnetic field magnitude at the location where the ion crystal 
is trapped above the surface. For the applications discussed below, one may require a 
magnetic field gradient along the trap axis, preferably high and constant over the few 
100 /im extent of an ion crystal. For the simulations [36 37 we assume a piecewise 
homogeneous current density according to the geometry of the center electrode (s). 
The results of our numerical simulations, presented in Fig. [5) show the magnetic 
gradient components plotted as function of z. The gradient component in z direction 
peaks with a width of about 100 /im at two locations at a distance given by the 
loop geometry. The transversal x-component of the gradient becomes constant, when 
moving the ion along the trap axis out of the chip center. For a current of Iwi = 4 A 
and = — 10 A one can achieve a much higher value, up to 23 T/m at a magnetic 
field strength about 7.6 mT. The ^/-component is even higher with more than 50 T/m. 

In the experiment, the magnetic field gradient will be switched on only for a short 
time such that the resistive heat load is reduced. In Fig. |6] we present current tests 
with the first generation of traps. The ratio of the measured voltage drop over the trap 
to the current determined via the voltage drop over a shunt is plotted as an increasing 
resistance over time. For a current of 6 Ampere, corresponding to a current density of 
4 • 10^ A/cm^, we are able to apply the load for more than 0.1 s. For larger currents, 
we observe a faster resistance increase since the resistance of the gold wire on the chip 
increases with temperature. It has been found that this ultimately leads to melting of 
the structure [38j. For transmitting such high currents, we use multiple bonds where 
up to 10 bonding wires of diameter 25 /im are used in parallel. Our tests show that 
only then the gold layer on the chip is destroyed before the bonding wires break. 
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Fig. 5. Magnetic gradient vl^l ^5 2/ ^ind z direction. The subplot a) shows the current 
flow within the two wires Wi and W2. b) The upper limit for the transversal gradient in 
our trap is reached with Iwi = 4 A and = — 10 A. c) Iwi — —Iw2 = 10 A results in a 
higher axial gradient, which is shown magnified in subplot d). 



In future, one might improve the cooling of the sapphire substrate, e.g. by a stack 
of Peltier elements, for even higher current densities: reducing the temperature by 
5T ^ 100 K would increase the gold conductivity by more than 30 %. It is reported 
that only an increase of more than 50 % in resistance will lead to permanent damage 
of the conductor [38 . For the planned experiments one would require a few ms spin- 
spin interaction times, such that more than 10 A appears within reach. A high current 
in the central conductor goes along with a certain voltage drop, which will have an 
impact on the z-position of the ion crystal if not counteracted by a suited voltage 
bias on the current supply. 
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2 Spin-spin interactions with ion crystals via the Stern-Gerlach 



Ions of mass m and charge q are trapped in the combined dc and ac electric potential 
resulting in secular motional frequencies cj/(27r) of typically a few MHz. With N ions 
bound in one crystal 3A^ modes of vibrations exist, and we will review this theory in 
the first part of this section. If the spin- 1/2 particles are exposed to a magnetic gra- 
dient, the position of the ions will change slightly due to a new force balance between 
the magnetic gradient-induced force on the spins and the electric restoring force of the 
harmonic trap potential. For a multi-ion crystal this results in an effective spin-spin 
interaction between the n^^ and the m^^ ion, denoted by Jn,m The weaker the har- 
monic restoring force, the more the ion positions are modified by the applied magnetic 
gradient and the larger are the mutual spin-spin interactions. We recapitulate briefly 
the interactions for a magnetic gradient along the 2;-axis which result in interaction 
strengths typically well below 1 kHz. Aiming for stronger couplings for quantum sim- 
ulations would require significantly stronger magnetic gradients. In this context we 
emphasize the usefulness of transversal magnetic gradients on linear crystals [21 39 . 
since we find that the spin-spin interactions largely increase when the trap control 
voltages are tuned such that the ion crystal is close to a structural transition insta- 
bility. We show that under these conditions, the corresponding spin-dependent forces 
even overcome the trap restoring force. The resulting interactions Jn^m are larger by 
orders of magnitude as compared with axial coupling strengths. In this case the ion 
crystal will rearrange, contradicting the initial assumptions of the calculation. The 
last part of this section outlines time dependent magnetic gradients with a current 
I{t) = /osin(a;t), where a resonant enhancement is employed to increase and design 
the couplings Jn,m- 

The total potential is the sum of the trap potential, approximated by a harmonic 
shape of the ions excursions Xn,i about their equilibrium positions in cartesian coor- 
dinates i and the electrostatic mutual Coulomb repulsions between the ions in the 
linear string: 



where = /AneQifnujI denotes a scaling factor and = (joI/ujI denotes the anisotropy 
between the axial and radial confinement. The Hessian matrix M (see Apendix [A|) 
determines all motional degrees of freedom: The eigenvalues of M are the squares 
of the motional frequencies 2TTVn of the ions, while eigenvectors of M indicate the 
directions of the ion motion. As long as one stays below a critical value of a^, ions 
align in a linear string and M reduces to a block diagonal matrix [40 41^ with en- 
tries muo1{B^ ^ By ^ A) . The eigenfrequencies of can be expressed by those of A: 
yJXjoii + 1/2 — (2ixv^Y' 12. As one can see, there is a critical value for the anisotropy, 
given through acrit = 2/((27ri^)^ — 1) ^ 2.53A^~^-^^ [42 , where v is the highest eigen- 
frequency of Matrix A. When crossing this point, the lowest eigenfrequency of B^ 
becomes imaginary and the linear chain is no longer stable. By increasing ol further 
the ions will rearrange in a zig-zag configuration to minimize the electrostatic energy. 



^ The MAGIC interaction employs Stern-Gerlach forces, but for a multi-ion crystal. Similar 
magnetic forces on a single electron or positron stored in a Penning trap were used for 
detecting flips of its spin orientation as in the seminal experiments by H. Dehmelt for the 
measurements of the g-factor 
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With an additional magnetic field gradient the Hamiltonian in one dimension for 
a linear string of N ions is given by [IE]: 



N 



N 



E 



Jn,m^z,n^z 



The first term describes the internal energy of a two level system due to the Zeeman 
shift hujn/'^ for the n^^ ion at its equilibrium position via the Pauli matrix a^^n- The 
confinement in the electric potential leads to vibrational eigenmodes with frequency 
Un- The last part represents a pairwise spin-spin coupling between the ions: 



h duJn 



2 dxji^i 



duJn 



determines the coupling strength using the Hessian M of the potential. The reason 
for such spin dependent forces is that the position of ion m is slightly modified when 
flipping an ion's spin, leading to changes in the Coulomb repulsion energy on ion n, 
and therefore changing its position also. Ion m experiences now a different Zeeman 
energy. Positive values of J result in an energy reduction for parallel spin orientations, 
whereas negative values lead to reduced energy for anti-parallel spins. Decoherence 
effects are expected from fluctuations of the magnetic fields and might be reduced 
by spin echo sequences as low frequency Fourier components typically dominate the 
noise-spectrum. 

J / 2n [Hz] 
0.15 



Fig. 7. Coupling constants Jn,m in a 
ten ion crystal at Uz = 27t ■ 310 kHz, 
for a b=l T/m magnetic field gra- 
dient in the z-direction, thus aligned 
along the ion crystal. A maximum 
Jn,m of 27V • 0.15 Hz/T is achieved. 




2.1 Magnetic gradient along the trap axis 

A magnetic field along the trap axis B = (60 + 2: • b)ez results in a field bo, which 
provides a quantization axis and which splits the \m = +1/2) and the |m = —1/2) 
spin states by an energy larger than the vibration quanta energy. The second term 
z ■ b results in a position dependent Zeeman shift, leading to forces against the axial 
vibrational modes only which are determined by the block A in the Hesse matrix M 
yielding: 

rax 9j^B ^ / 4-lx /rx 

with Bohr magneton \ib and Lande factor gj = 2. The relative magnitude between 
different entries of J depends only on the number of ions, whereas the absolute height 
is proportional to b!^ /oo^. Fig. [t] shows the coupling constants for ten ions, normalized 
to a unit magnetic field gradient of 1 T/m . All couplings are positive and strongest 
for nearest neighbors with a strength of 2tt • 7.4 Hz for 7 T/m and an axial frequency 
of uj^ = 271- 310 kHz. 
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2.2 Magnetic gradient transverse to the trap axis 

We now investigate the interaction when the gradient is aligned along the x-axis and 
the spin dependent force acts against the transversal eigenmodes which are governed 
by the anisotropy parameter. The transverse confinement and the lowest radial eigen- 
frequency can be reduced without modifying the distances between the ions. The 
Stern- Gerlach force acts by moving the ions transversally away from the trap axis, 
and the lowest frequency transversal eigenmode, the zig-zag mode governs the inter- 
action. This mode is excited most efficiently with an alternating, anti-ferromagnetic 
spin alignment. The transversal magnetic gradient couplings Jn,m are determined by 
the matrix (see Appendix [A| and shown in Fig. [s] for different values of a. When a 
approaches its critical value, we find coupling strengths of about 27r • 158 kHz, where 
the frequency of the zig-zag mode Uzz is at 11 kHz and the transversal field gradient 
is 23 T/m. Whereas for values far from the zig-zag instability the pattern displays 
mostly next neighbor negative couplings, the interactions increase when approaching 
the instability and show a fully antiferromagnetic behavior. The strength of mutual 
interactions is stronger in the middle of the ion crystal, in contrast to the behavior 
of the Jm,n in the case of axial magnetic gradients in Fig. [Tj The radial restoring 
forces in an ion crystal are weakest in the middle, while the axial restoring force is 
weakest at the ends. We stress that the coupling strengths with transversal magnetic 
gradients are about three orders of magnitude higher than in the case of interactions 
driven by a magnetic field gradient of the same magnitude but in the z direction, 
when lowest radial eigenmode frequency is reduced to about 10 kHz. Gradients in 
the radial direction are generated by the current running in the central wire Wi (see 



5). An alternative scheme for realizing a transversal magnetic field is sketched in 



Fig. 

Fig. |9j A second chip beneath the trap containing just two broad wires with counter 
propagating currents (narrow U-shaped structure), creates a vertical magnetic field. 
Working with the same current density as before will give rise to a gradient of about 
40 T/m at a height of 400 /im over the plane of the U, fitting perfectly to the ions 
location. The axial component vanishes, as the bent part of the current loop is far 
from the ion position. With the above values we calculate a coupling strength of up to 
27T ■ 479 kHz for u^z = 11 kHz. One question to be answered by future experiments is 
how accurately the trap control parameters can be adjusted in order to approach the 
instability point even further for increasing the couplings even into the MHz regime. 
The MAGIC interaction does not strictly require cooling to the vibrational ground 
state. Still, at small values of Uzz the mean phonon number at the Doppler tempera- 
ture will be excessively high. One possible way to reduce the phonon number in the 
mode, that mostly contributes to spin-spin coupling, would be to cool that mode still 
at high Uzz and then adiabatically ramp up the value of a for achieving the desired 
spin-spin coupling strengths. 



2.3 Initialization and read-out of ion spins 

Before the application of the spin-dependent force, one might like to align all spins 
in one direction. This can be done by standard optical pumping e.g. with (j+ laser 
light such that the m=+l/2= | t) state is occupied. An even higher fidelity of the 
initialization employs frequency-resolved optical pumping yielding a transfer of 99.6% 
[43] . A subsequent optical Raman pulse on the entire ion crystal with a well chosen 
length of 6 = 7r/2 generates the superposition state (| t) + I i))/V^ for each ion. 
The read-out would be performed by a detection of state-dependent fluorescence, a 
well established tool in experimental quantum information with trapped ions. The 
advantage is that, with careful application, the read-out fidelity can reach values 
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C) 1 2 3 4 5 6 7 8 9 10 J / 2n [Hz] 




Fig. 8. Coupling constant of ten ions normalized to b=l T/m magnetic field gradient in 
the radial direction when approaching the critical anisotropy for ten ions (acrit = 0.047348). 
The three subplots indicate in color coding the transition into a ferromagnetic coupling with 
red for positive and blue for negative values of J. a) Choosing a value of a = 0.0097819 
results in z/^^=2.8 MHz and a maximum Jn,m of 27r • 7.1 10~^llz/T. b) For a = 0.032669, 
jy^^ = 1.7 MHz and Jn,m increases to 27r • 7.1 10~^Hz/T. c) With a = 0.047347 we find 
=10.7 kHz and Jn,m of 27r-300 Hz/T. 




Fig. 9. Sketch of a planar micro ion trap for 
applying radial magnetic forces: a) Trap chip 
for ac and dc electric potentials, b) underly- 
ing second chip to be mounted directly under- 
neath (a) for magnetic field gradients in the 
radial direction. Arrows indicate how (b) is 



mounted under (a). 
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exceeding 99.99% [47], and that this technique also resolves the individual ion sites 
such that the spin order after the interaction can be directly observed by state- 
dependent fluorescence on CCD-camera images. 

2.4 Selective excitations via resonant current 

So far, we have discussed spin-spin interactions when a static magnetic field gradient 
is applied. The small induction of the wire loop that generates the gradient allows one 
to apply alternating currents I{t) = /osin(a;t), with frequencies of a few kHz up to 
even a few MHz and where uj may be tuned close to one of the motional frequencies 
27ri/^ in the ion crystal. After a proper initialization, this will predominantly excite 
the ions at the positions corresponding with the eigenvectors of that specific mode 
m, thus selecting spins at specific positions of the crystal |44|45j . If for example, the 
breathing mode in a 3-ion crystal at {27r)V3uz is excited by an alternating gradient 
in the z-direction, only the outermost ions couple, and the middle ion is excluded 
from the interaction. When multiple frequencies are applied on ion crystals with 
I{t) = Ij sm{ujjt)^ one may excite superpositions of eigenmodes, where arbitrary 
ion excitations can be achieved. 

The second advantage of resonantly applied gradients results from a resonant en- 
hancement, such that the interaction strengths Jn,m are increased as compared to 
their static values. We propose the generation of Schrodinger cat states [46] \^f) = 
^{\ t, + i| I, — where the ion crystal is first initialized into the superposi- 

tion spin state and, via the magnetic gradient forces, the motional degree of freedom 
is entangled with the spin state of an ion when a single ion is excited very close to 
the axial vibration mode with a detuning 5 = uj — 27rUz. Large values of displace- 
ment a are reached for small S. In analogy to the application of spin dependent light 
forces a two-ion crystal gate operation similar to the geometric gate operation with 
state dependent ac light forces [48 l49|5Q|51|52j can be realized [5T, where a transient 
motional excitation is generated only for specific spin (qubit) states, and results in a 
conditional phase shift of the wave function. 

It should be noted that the investigation of the linear to zig-zag configuration 
in a linear string |54l55j is a quantum phase transition which has attracted much 
attention recently. However, the experimental realization has not yet been achieved. 
We state that the transverse magnetic field gradients would be ideally suited to drive 
the system in a spin-dependent way from one phase to the other [56] . 




Fig. 10. Sketch of a planar trap for trapping 
multiple parallel ion strings. 



3 Spin-Spin interactions in two dimensional ion crystals 



Having discussed the tunability of spin interactions with linear ion crystals, we aim for 
simulations of larger systems of interacting spins and more complex geometries, that 
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are only possible in two dimensional systems. In the following, we will outline in detail, 
how planar ion traps with integrated magnetic field gradients can reach such advanced 
goals of quantum simulation. For this we present a calculation of vibrational modes 
for 2D-ion crystals, followed by the results of spin-spin interactions in such system. 
Our versatile, dynamic and even possibly fast control of the ion crystals positions 
in the two-dimensional lattice, can be seen in contrast to proposals for static planar 
arrays [23 57 . It appears that the advantage of our scheme is that it may allow for 
investigating frustrated spin systems with tunable interactions. The geometry of the 
latices can be changed on the run. This may be also reached by a dynamical control of 
the radio frequency, as proposed recently ^58^. Furthermore, our scheme gives a clear 
route how the ions crystals are loaded in a remote, unproblematic since wider part 
of the trap and then being transported to the narrow trap interaction region for the 
two dimensional spin-spin interactions. As the loading of linear shaped planar traps 
has been successful in many labs including ours, this scheme appears to be in reach 
of the experimental possibilities. 

A sketch of a planar trap for the storage of two linear ion strings parallel to 
each other is shown in Fig. [lOj extending the generic planar five-wire trap design 
only slightly, here by one rf and one dc electrode. Lithographically produced micro 
structures [59 will allow for trapping multiple strings. For describing the interactions, 
the potential in eq. |4]has to be slightly modified by adding the electrostatic coupling 
between ions which sit in different strings, 

where N now denotes the number of ions in one chain and d the distance between 
the trap centers in the x-direction. The upper limit of the sums in the second part of 
eq. |4] therefore runs to 2N . The general form of the Hessian M changes to dimension 
2A/', and M is still blocked, but no longer diagonal, for there are non-zero entries 
coupling radial and axial modes. Only in the limit of d to infinity does M simplify 
to a diagonal matrix with separated entries {B^ , ^ ^ ^ A) describing two 
independently oscillating linear strings. 



3.1 Eigenmodes of parallel arranged linear ion crystals 

In a two-dimensional ion crystal, any movement of one chain will influence the other 
such that the separation in radial and axial eigenmodes fails, and 3x(2A^) eigenmodes 
exist. In the following we exemplify the effects with a system of two strings each con- 
taining two ions. Ions in one string can vibrate in common or contrariwise, commonly 
denoted as center-of-mass and breathing mode respectively. The respective frequen- 
cies u are characterized by the square root of the Hessian matrix's eigenvalues given 
in multiples of mcj^, namely 1 and a/3 in the axial direction. For the radial direction 
mode frequencies depend on the value of a, and are denoted as the radial center-of- 
mass and the rocking mode. In combination with the second string arranged parallel 
to the first one, we find that the four modes of one string split up in a co- and an 
anti-mode, plotted over the distance d in Fig. [IT] (tables of eigenmodes in the app.[T]). 
The frequency shifts are stronger if the chains are moving in the anti-configurations. 
The modes with dominant amplitude along the x-axis in the direction of d display 
increasing eigenvalues, whereas the z dominated transversal modes are decreasing in 
frequency. With many ions per chain and below a certain a both chains will undergo, 
similar to the case of a single string, a transition to a zig-zag-zig-zag configuration. 
However at small values of d both chains generate a stabilizing effect onto each other 
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Fig. 11. Eigenfrequencies v of 
two ion strings with two ions 
each confined in a 2-dimensional 
array, plotted in multiples of 
uJz/{2ti) as a function of the dis- 
tance d between the axes of the 
strings for N — 2. The values 
correspond to a radial confine- 
ment with 0Jx,y,z of 3.42 MHz, 
2.72 MHz, and 1.04 MHz, when 
d approaches infinity. The plot 
shows the dependence of both 
radial and both rocking modes 
(green), the breathing (red) and 
the axial center-of-mass mode 
(green), from top to bottom. 
The anti-modes are plotted in 
bright color and the co-modes in 
dark, respectively. The ion's mo- 
tion in one string is indicated by 
the direction of the arrows for 
large distances. 
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due to the mutual Coulomb repulsion, thus remaining in a linear configuration well 
beyond the critical value of a for a single ion crystal of same size. We also calcu- 
lated the resulting equilibrium positions for the two-dimensional crystal and observe 
a small bending of the strings due to a comparatively greater shifting of the ions in 
the middle of the strings. In the case of ten ions at a distance of d = 50 jam with 
ax = 0.00975, the numerical solution yields a relative increase of the distance between 
the strings of 1.4 10~^ and a relative increase of the inter- ion distances in one string 
of 1.4 10~^. This effect might be counteracted by proper compensation. 



3.2 Magnetic coupling with axial and radial gradient fields 

Surprisingly complex magnetic interactions arise from the application of magnetic 
field gradients. In Fig. [12] the interaction strengths are shown for two chains with ten 
ions with a normalized 1 T/m magnetic field gradient applied either purely in z or 
purely in the x direction. The coupling matrix is separated into four blocks, where the 
diagonal ones represent the intra chain coupling, as discussed in the previous sections. 
The off-diagonal blocks describe the interaction between chains. 

In the case of a field gradient aligned along the direction of the ion strings (in 



Fig. 12 1), the ferromagnetic coupling within a string is recovered. A very weak cou- 
pling to the neighbor string is anti- ferromagnetic, and thus the lowest energy state of 
the system will lead to two spin chains, where all spins within a string are parallel, but 
oppositely oriented with respect to the other string. For b = 40 T/m the maximum 
magnitude of J is 27r-260 Hz within one string and 27r-7.2 Hz between strings. Even 
more interesting behavior is observed in the cases displayed in Fig. [12] b,c,d for the 
transverse gradient. When approaching the zig-zag instability, the coupling within 
the strings develops an anti-ferromagnetic order, while the coupling between different 
strings undergoes a transition from ferromagnetic to anti-ferromagnetic order. Far 
above the critical value of a, at a = 0.009782 (in Fig. [T2]3) we calculate very small 
values of J with 27r-l.l Hz and 27r- 2.4-10"^ Hz. Approaching closer to the instability 
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Fig. 12. Coupling strength for two chains with ten ions at a distance of 50 /xm normalized to 
b = 1 T/m unit magnetic gradient. The inter-string coupling blocks (off diagonal elements) 
are found in a 3D plot below the corresponding matrix plot. The gradient magnetic field 
is applied either in the z (axial) direction in case of subplot (a) or in x (transversal) for 
the subplots (b,c,d). In (a) the anisotropy parameter is fixed well above the instability to 
the zig-zag transition at a = 0.009782 with uoz = 310 kHz yielding ferromagnetic coupling 
within the strings and much weaker anti- ferromagnetic interactions between both strings. 
For a transverse magnetic gradient, the anisotropy parameter is varied from a = 0.009782, 
Uzz = 2.8 MHz (in b) to a = 0.047347, a value which corresponds to the instability transition 
for a single linear string, at i^zz—^^Q kHz (in c) and a — 0.048014 with Vzz = 11.5 kHz (in 

d). 
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we find 27r-416 Hz and 27r • 2.9 • 10~^ Hz (in c) and when the zig-zag mode is tuned 
down to = 11.5 kHz the couplings increase to 27r-402 kHz and 27r • 2.1 • 10~^ Hz. 
Even though the off-diagonal elements of the coupling matrix are still rather weak 
as compared with the diagonal elements, we expect that the technique of alternating 
near-resonant field gradients as discussed in Sec. |2.4| will enhance these elements such 
that an inter-string coupling is possible to observe. Evidently, bringing the two ion 
strings to a distance even closer than (i = 50 /im would also lead to significantly 
increased off-diagonal couplings Jn,m- 



3.3 Frustration effects 



A configuration of two strings parallel to each other may also be used to investigate 
frustration effects in a spin system. For this we propose shifting one of the chains 
along the z axis by half the inter-ion distance in the string. Now, ions of one string 
are at equal distance to two ions in the other chain. With a ferromagnetic interaction 
between both strings, and an anti- ferromagnetic interaction within the strings, the 
ions minimize energy when aligning with two counterparts in the other string, but 
those have different spins, as shown in Fig. 13 If the bending of the strings can't 
be overcome, this would be a major drawback to this scheme, since in that case 
only one triple of ions can be exactly frustrated at one time. The calculation of 
the resulting spin orientation is beyond the scope of this paper, but in principle it 
is possible to adjust the intra- and inter-string coupling strengths. As an example, 
by bringing the strings as close as 20 /im together and with a relatively loose axial 
confinement of cj^ = 27r • 155 kHz and a = 0.0024 the intra-chain interaction strength 
drops to 27T • 1.02 • 10~^ Hz, whereas the inter-chain coupling strengths are about 
27r • 2.83 • 10~^ Hz, thus leading to a tunable ratio, in this case of about 3.6. In a 
future experiment, one may control the mutual interactions by a variation of a and 
the magnetic field gradient, also by positioning of the strings with small changes 
of the trap control voltages such that it appears well within reach to switch back 
and forth between rectangular lattice configuration and a triangular one resulting in 
frustration effects. The ion model system might allow one to contribute to the vivid 
theoretical discussion of quasi one dimensional structures in interacting quantum 
wires [60 , where the quantum conductance depends on the filling factor, as when 
increasing the number of electrons in a quantum wire, they undergo a transition to a 
zig-zag configuration. It may be further investigated how trapped ion crystals could 
be employed to simulate transport in quasi-linear solid state structures. One might 
use trapped ion strings which are moved along each other as a model system with 
versatile and tunable properties, as a model of genuine solid state physics problems. 
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Fig. 13. Displacement of two linear ion strings 
with respect to each other creates a triangu- 
lar geometry. The designed coupling of spins 
could be used to create frustrated ion crystals. 
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4 Conclusion 

We have outlined a scheme for investigating spin-spin interaction with the help of 
planar ion traps. We obtained a description of the magnetic interaction properties for 
a one or two dimensional crystal, strongly depending on both the direction and the 
strength of the magnetic gradient that is applied. Ferromagnetic or anti-ferromagnetic 
configurations are generated, and with a control of the gradient by two independent 
current-carrying wires as sources for the fields it is possible to direct it in any direction 
of the x-z plane. The strength of the interaction is proportional to two easily accessible 
parameters, namely the magnetic gradient via the chosen current and the axial trap 
frequency determined by the applied dc voltages. Although the electrode geometry of 
a planar trap is fixed, it still offers a range of possible trapping frequencies to realize 
potential wells that confine the ions very tightly or allowing the crystal to restructure 
itself e.g. into the zig-zag configuration. The segmented outer dc electrodes can be 
used to displace ion strings against each other or to split up one chain. Offsets on the 
rf electrodes could be used to modify the distance between chains, thus one can shift 
the eigenfrequencies of the system and manipulate the inter chain coupling. A rich 
plethora of interactions is possible in such two-dimensional quantum spin crystals and 
the related physics opens up a wide field of multi-particle spin simulations. Perhaps 
among the most striking features is the possibility to simulate frustrated spins in 
two- dimensions. In the paper we have discussed all these future possibilities in very 
close correspondence to the existing experimental parameters, thus in a fully realistic 
setting. 
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A Calculation of the Hesse Matrix 



Since we have N ions and in general three spatial degrees of freedom the Hessian of 
Eq. [4] can be written formally as 



jN-\-o,kN-\-p 



k ^op 



an 



3(^71 



2 ^ 



(So 



^mp^ 



with o and p denoting the ion's number, j and k the spatial direction and 6jk the 
Kronecker delta. For a single string and as long as one stays below the critical value 
of a, ions align linearly and thus M reduces to a block diagonal matrix [41p40) with 
entries muj'^{B^ , ^ A) where 



1 + 2^^=1^;;:^ ifn = m 
I- _- 13 if n 7^ m 

I Zm Zn I 



In the case of two strings with each holding N ions, the two sums in M can be 
split apart in intra string and inter string components: 



2N 2N 



N 



2N 



N 



2N 



EE = E + E +25: ^ = (a) + (6) + (c). 



n=l m=l n,m=l 



n,m=A/'+l 



n=l m=N+l 



Sum (a) will yield to the matrices A and B for one chain, (6) to that of the other 
chain, such that one has a matrix (5^^, 5^^, 5^^, 5^^, A^ ^A^). Part (c) however leads 
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to new entries in M, which reads with = {Sno — ^mo){^np — ^mp)- 



op ~ ^ / J ^nm 



-nm Ix -X |5 
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^22; _ 73 rop 



™ Ix -X |5 

n,m 

/-^xz _ 73 rop ^{^m ~ ^n)(^n ~ ^m) _ ^zx 
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B Trapping parameters and eigenmodes 

Table ([T]) shows possible values for the voltages applied to three neighboring pairs of 
dc electrodes to achieve axial confinement. The electrode pair in the middle is held 
at voltages U2 and both outer pairs at [/i. 



Table 1. Influence of the axial confinement on radial frequencies and resulting anisotropy. 
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a;^/27r[MHz] 
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cj^/27r[MHz] 




ay 


1.30 
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-6.13 
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0.0542 


0.0446 


13.0 


-12.2 


2.77 


3.39 


0.98 


0.125 


0.0835 


23.2 


-21.7 


2.41 


3.55 


1.30 


0.297 


0.135 



The following tables contain the eigenmodes of a system with two ion chains, each 
containing 2 ions. The subscript refers to the nth ion, where ion 1 and 2 are in one 
chain and ion 3 and 4 are in the second chain. The different modes are sorted in 
descending order according to the eigenvalues. 
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Table 2. Eigenmodes in a distance d = 100 //m. 
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Table 3. Eigenmodes in 
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Table 4. Eigenmodes in a distance d = 20 /xm. 
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